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استفاده از روش‌های نظریه گراف‌ها برای بهبود خوش وضعی ماتربس نرمی سازه ها * 
فرزاد شفیعر و ۱ مصطف خانزادی ٩(‏ 


چکیده در تحلیل رایانه‌ای سازه‌هاء تحلیل باید به‌صورتی انجام پذیرد که علاوه ب رکاهش فضای کامپیوتر و زمان اجرای محاسباتی کامپیوتر 
عطاهای محاسباتی نیز به کم‌ترین مقدار ممکن برسد. این امر با تشکیل ماتریس‌های سازه‌ای مناسب که در حین پر صفر و وش ساختار بودن» 
حوش حالت نیز می‌باشد, تحفق خواهد یافت. این یک مسأله‌ی بهینه‌سازی, چندمنظوره خواهد بود. در اين مقاله بر روی خعطاهای محاسباتی, اعم 
از عطاهای گرد کردن و ... در تشکیل ماتریس نرمی سازه و کاهش این عطاها با ارائه الگوریتم‌های مناسب, پایدار و کارا تأکید شده است. زمانی 
که به دلیل ساعتارهای به‌نحصوصی, ماتریس نرمی به‌صورت بد وضع درم یآید عطاهای محاسبات افزايش می‌يابند. این حالت در سازه‌های دارای 
اعضای با سخنی‌های متفاوت مشهودتر است. استفاده از اعضای ضعیف, سبب ایجاد درایه‌های بزرگ غیر قطری در ماتریس نرمی می‌شود که 
نتیجه ی آن ایجاد بد وضعی ماتریس نرمی است. در محاسبات عددی, بد وضعی یک ماتریس را با جابه‌جایی سطرها و یا جایگزین کردن چند 
میم یوت یراق نان توف یاب ام اشت :که تراسا خری ما از تشکین مالوینی ترمی دصر قاه لاعفا کرت رای 
این منظور از نحصوصیات نوپولوژیکی و جبری ساده در تحلیل, به طور موثر استفاده می‌گردد و این امر با معرفی پایه چرخحه‌هایی با حواص وه 
صورت می‌پذیرد که منجر به تشکیل ماتریس نرمی با شرایط موزونی شبه بهینه شده و در نتيجه عطاهای محاسباتی به مقدار قابل توجه ی کاهش 


می‌یابد. 


واژه‌های کلیدی نظریه گراف‌ها. روش نیرو. خوش وضعی ماتریس‌های سازه‌ای 
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۳ نسخه‌ی اول مقاله در تاریخ ۸۱/۷/۲۷ و نسخه‌ی نهایی آن در تاریخ ۸۸/۹/۸ به دفتر نشریه رسیده است. 
(۱) نویسنده‌ی مسئول. کارشناسی ارشد. دانشکده‌ی مهندسی عمران, دانشگاه علم و صنعت ایران, نارمک. تهران 


(۲) استادیا دانشکده‌ی مهندسی عمران دانشگاه علم و صنعت ایران, نارمک. تهران 
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۳/۸ 


مقدمه 

استفاده از رایانه‌های عددی برای حل مسائل سازه‌ای نیاز 
به حل یک سیستم بزرگ معادله‌ی جبری به‌صورت زیر 
است: 
)۱ 0 < :۸ 

معادله‌ی (۱) هم برای روش نیروها و هم برای 
روش تغییر مکان‌ها باید حل گردد. بعضی وقت‌ها جواب 
معادله‌ی (۱) با تغییرات اندک در درایه‌های ماتریس ۸ 
ویا 0 ویا هرد تغییرقابل توجهی می‌کند. در این صورت 
گفته می‌شود ماتریس ۸ نسبت به این حل. بد حالت 
است. دقت حل معادله‌ی (۱) می‌تواند به مشخصات 
اون سا شسپیای فان مرها نش ان 
مشخصات و رابطه‌ی آن‌ها با منبع انتشار و پخش 
خطاهای ممکن حائز اهمیت است. بدین منظور 
روش‌های بهتر فرمول‌بندی مسأله باید پیدا شود و 
تکنیک‌هایی برای پیشگویی.کشف وکمینه کردن خطاهای 
حل, باید مورد توجه قرار گیرد. بد حالتی ماتریس نرمی 
و ماتریس سختی سازه برای روش نیروها و روش سختی 
توسط کاوه [1-2] مطالعه شده است. هم‌چنین بد حالتی 
ماتریس سختی سازه برای روش سختی به وسیله‌ی شاه 
[3] مطالعه شده است. در کار وی روش‌هایی برای بهبود 
خوش وضعی ماتریس سختی سازه داده شده است. یک 
تحقیق ریاضی از تحلیل خطای ماتریس مربوط به 
وتان وگن که | است: بر کار انفان فان واه 
کم کا اد لاش تابابد ای یکین اتسار سا 
پایدار فیزیکی ناشی شوند. اثر زیر سازه‌ها روی شرایط 
ماتریس سختی به‌وسیله‌ی گرومز و راو [5] تحقیق شده 
این نتیجه به‌دست آمده که زیر سازه‌ها به طور عمده روی 
دقت حل سیستم‌های بد حالت تأثیری ندارند. فیله و [6] 
یک روش متعامد سازی برای بهترین شرایط ماتریس 
نرمی وسختی پيشنهاد کرد؛ هر چند این یک روش غیر 
عملی برای سازه‌های پیچیده با اعضای زیاد است. بهینه 
سازی شرایط معادلات تعادل, درزمانی‌ که یک روش 


جبری به‌کار رود. به‌وسیله‌ی رابینسون و هاگنماخر [7] 
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استفاده از روش‌های نظریه گراف‌ها برای بهبود ... 


مطالعه شده است. اخیراً روش‌هایی برای انتخاب انواع 
ویژه پایه‌های استاتیکی و سینماتیکی که منجر به خوش 
حالتی ماتریس نرمی و سختی سازه نسبت به پایه‌های 
کلاسیک گردیده است. توسط کاوه 1 توش ی یافته 


استتا: 


روش نیروها 
تحلیل سازه‌ها به روش نیروها که در آن نیروهای اعضاء 
به‌عنوان مجهول به‌کار می‌روند. از این نظر جذاب است 
که خواص اعضاء یک سازه و طراحی آن‌ها غالبا به 
نیروهای اعضاء بیشتر از جابه‌جایی گره‌ها بستگی دارد. 
انم روش اسان ۱۹ ره وی داهتمبیا تور 
رایانه و مناسب بودن روش تغییر مکان‌ها برای رایانه‌ای 
شدن. اکثر محققین مجذوب روش تغییر مکان‌ها شدند و 
در نتیجه روش نیرو و بهینه کردن آن دچار فراموشی شد. 
چهار طریق مختلف برای تحلیل سازه‌هابه روش 
روا فان شل اشست که ها یل ۱ 
روش‌های توپولوژیکی. 
۲ روش‌های جبری. 
۴ روش‌های مختلط توپولوژیکی- جبری» 
12 روش مجتمع. 
روش‌های توپولوژیکی به‌وسیله‌ی مان در [9] 
هندرسن [8] برای اسکلت‌های ساختمانی با گره صلب با 
به‌ کار بردن انتخاب دستی برای پایه چرخه‌های مدل 
گراف آن‌ها توسعه داده شده است. روش‌های مناسبی 
برای رایانه‌ای کردن آن‌ها توسط کاوه [10,11] ارائه شده 
است. این روش‌ه ای توپولوژیکی برای کلیهی 
اسکلت‌های ساختمانی تعمیم داده شده است. به‌عنوان 
مثال» برای قاب‌های صلب. خرپاهای مستوی مفصلی و 
خرپاهای فضایی [11]. 
روش جبری توسط کانکو وهم‌کاران [15] تویسو 
[14] رابینسون وهم‌کارانش [13] دنکه [12] توسعه داده 
شده است و روش‌های مختلط توپولوژیکی- جبری 
به‌وسیله‌ی گیلبرت و هم‌کارانش[16]؛ کلمن و پوتن 


سال بیست و دوم شماره یک» ۱۳۸۹ 


فرزاد شفیعی دیزجی- مصطفی حانزادی 


قفا ورین [ تلا بهکان کدف شاه اس 
روش مجتمع نیروها به‌وسیله‌ی پتنایک [20-22] 
توسعه داده شده است که در آن معادلات تعادل و شرایط 
سازگاری به طور هم‌زمان برحسب متغیرهای نیرویی بیان 
شده‌اند. در این فصل. روش‌های توپولوژیکی که از 
مفاهیم کلی تثوری گراف‌ها استفاده می‌کنند اراشه شده 
اسر ناف یت فزرشوال یی کته پر ناش 
معادلات تعادل نیروهاء سازگاری تغییر مکان‌ها و روابط 
بین نیروها و تغییر شکل‌ها استوار است؛ شرح داده 
می‌شود. یک سازه ۸ با عضو ۷ و لا گره که (۷65 
درجه‌ی نامعینی استاتیکی ([5) دارد را در نظر بگیرید. 
تعداد 76٩(‏ نیروی مجهول مستقل را به‌عنوان مجهول 
اضافی انتخاب کنید. این نیروهای مجهول را می‌توان از 
عکس‌العمل‌های خارجی و یا نیروهای داخلی سازه 
انتخاب کرد. مجهول‌های اضافی به‌صورت زیر تعریف 
می‌شوند: 
(۲( 5 ۹ 
قیود مربوط به مجهول‌های اضافی را برای تعیین 
یک سازه معین استاتیکی که به‌عنوان سازه پایه 5 (آزاد 
شده یا مقدماتی) شناخته می‌شود. از بین ببربد. بدیهی 
است که یک سازه پایه. بایستی پایدار باشد. بارهای 
گره‌ای را می‌توان با عبارت زیر تعریف کرد: 
)۳( تنج اک 
بدین علت که بارها فقط در بعضی از گره‌ها ممکن 
است اثر کنند» به‌جای  "(‏ به‌کار رفته است. حال نتایج 
توزیع تنش به سبب وارد آمدن بار 9 برای یک تحلیل 
خطی به‌وسیله‌ی روش یرو را می‌توان به‌صورت زیر 
)( +3 < ۲ 
که 32 و رظ ماتریس‌های مستطیلی شکل‌اند و هر 
کدام ۷ سطر و به ترتیب لا و (۷)5ستون دارند.۳۲ که 
به‌عنوان یک حل خصوصی شناخته می‌شود. تعادل را 
ارضا می‌کند و توسط بارها وضع می‌شود و 3:0 که حل 
مکمل یک مجموعه حداکثر از سیستم‌های خود متعادل 


سال بیست و دوم شماره یک. ۱۳۸۹ 


۳۹ 


تنش (5.۳.56) است» به‌عنوان یک پایه استاتیکی شناخته 
می‌شود. ستون‌های ۳ یک پایه‌ی استاتیکی ٩‏ را تشکیل 
می‌دهند. توجه کنید که بر روی هر چرخه سه سیستم 
خود متعادل برای قاب مستوی می‌توان تشکیل داد. 

یک سازه‌ی پایه لازم نیست به‌صورت یک سازه‌ی 
معین انتخاب شود. برای یک مجهول اضافی سازه پایه 
می‌توان اطلاعات لازم را يا به‌وسیله‌ی تحلیل اولیه برای 
بارهای ۲ و واکش‌های ((2,..,۵8با< )2 و یا 
به‌وسیله‌ی استفاده از اطلاعات موجود تهیه کرد. 

معادلات سازگاری در سازه‌های حقیقی حال 
ساخته می‌شود. با به‌کار بردن نسبت بار, تغییر مکان برای 
هر عضو و جمع‌آوری آن‌ها در قطر ماتریس نرمی مجزا 
ی 
(۵ 0رظر + م,رظیط ح تساحاا 

معادله‌ی بالا را می‌توان در یک شکل تقسیم‌بندی 


که اه 


۰ له سل 


که و۷ شامل تغییر مکان مربوط به اجزاء نیروهای ۲ 
می‌باشد و .۷ تغییر مکان نسبی برش‌ها را برای سازه پایه 


نم ۳ ,3,5 یط 
فا صصظ رظرظ ۲ 

یا 

)۸( 9( ,بوظ) + م(مظی۳نظ) ع ولا 

و 

)8( (ررترظ) + م(مظرظ) ۷ 


حال با اعمال شرایط سازگاری داریم: 


(۱۰) 20 ۷ 
با تلفیق معادله‌ی بالا با معادله‌ی )٩(‏ داریم: 
)۱( «(وظ,۳ظ) ‏ (رظرطظ)- و 


با جای گذاری در معادله‌ی (۸) داریم: 
)۱۲ رظ یط - رظر 3 ]ع و 


و نتایج تتش در یک سازه را می‌توان به‌صورت زیر 
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نوشت: 

۱۳( و[ رظ,۳ظ ۱ (ر۳,۳ظ) .ظ- رظ ]- ۲ 
که رظ ,۴ 13 < 0 به ماتریس نرمی یک سازه موسوم 

است. شرایط یک تحلیل بهینه به روش نیروها این است 

که ماتریس نرمی ) دارای شرایط زیر باشد: 

کت ماشنه 

۲ خوش حالت باشد. 

۴ خوش ساختار باشد (به‌عنوان مثال پهنای باند کمی 

داشته باشد). 

برای دسترسی به شرایط ۱ و ۲ ساختار ,ظ بایستی 
به خوبی طراحی شود؛ زیرا ساختار »۳ برای هر مسأله 
خاص قابل تغییر نیست. شرط سوم دارای ذات ترکیباتی 
می‌باشد. در این مقاله. خاصیت دوم مطالعه می‌شود. 

ماگزشی و که شام یایدهلی اتکی بی بافتته 
دارای الگوتی مشابه با ماتریس ) است (که ) ماتریس 
مجاورت چرخه عضو می‌باشد) به طور مشابه ماتریس 
,3 ,۲ 13 < 6 معادل بهینه کردن ماتریس همسایگی 
پایه‌های چرخه‌های گراف» 0 می‌باشد. این تبدیلات» 
مزایای بسیاری دارند از قبیل: 

مرتبه‌ی ماتریس 0 کسوچکتر از مرتبسه‌ی 6 
می‌باشد. برای یک سازه‌ی فضایی مرتبه‌ی 
ماتریس 00 شش بار. و برای یک سازه‌ی دوبعدی. 
سه بار کوچکتر از مرتبه‌ی ت) می‌باشد. 

۲ درایه‌هنای ساتریس‌هنای 6و 66 اعضای و7 
(مجموعه اعداد صحیح تشکیل شده از صفر ویک) 
می‌باشند. پس کار بر روی آن‌ها بسیار راحت‌تر از 
و ت) خواهد بود؛ زیرا درایه‌های ۳ و ت) اعداد 
تجهیتی هی بات 

۳ پیشرفت‌های به‌دست آمده در ریاضیات ترکیباتی و 
تئوری گراف به طور مستقیم قابل استفاده خواهند 
من 

۶ تناظر بین روش‌های تحلیل جبری و ترکیباتی 
پاپدارتر می‌ننود: 

همان‌طور که مشاهده شد برای دستیابی به یک 
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استفاده از روش‌های نظریه گراف‌ها برای بهبود ... 


تحلیل بهینه, مسأله‌ی تأمین شرایط تحلیل بهینه برای 6 
به مسأله تأمین شرایط مذکور برای 26 تبدیل می‌شود. 
این بدان معنی است که حتی قبل از این که به تحلیل یک 
سازه پرداخته شود با بررسی چرخه‌های پایه‌ی انتخاب 
شده می‌توان بهینگی تحلیل را مشخص نمود. از طرفی 
برای دستیابی به یک تحلیل بهینه» بایستی پایه‌ی 
چرخه‌های مناسبی انتخاب نمود. هم الگویی پایه 
استاتیکی یک سازه اسکلتی و یک پایه چرخه‌ی مدل 
گرافی آن دیده شد. روش‌های موثر برای انتخاب یک پایه 
چرخه هم الگو برای پرصفر کردن ماتریس 3 توسط 
کاوه [11-25-26]» کسل وهم‌کاران [24] هوبیکا و سیسلو 
[23] کنتریی یافیه است: الکوریت این بترای مرب 
کردن المان‌های پایه‌های انتخاب شده که منجر به کم 
عرضی ماتریس نرمی می‌شود توسط کاوه [27-28] 
مطالعه شده است. 


پایه چرخه‌های کمینه و بهینه 
همان‌طور که می‌دانیم. ماتریسی تنی نامیده می‌شود که 
تعداد زیادی از درایه‌های آن صفر باشد. اهمیت تنکی از 
این لحاظ است که کار کردن با ماتریس‌های تنک. باعث 
صرفه‌جویی در زمان محاسباتی رایانه می‌شود. از طرفی 
اکثر ماتریس‌هایی که در تحلیل سازه‌های عملی به‌دست 
می‌آیند می‌توانند تنک شوند؛ حتی اگر چنین نباشند. 
«ضریب تنکی» 7 برای یک ماتریس به‌صورت تعداد 


درایه‌های غیر صفر آن تعریف می‌شود. 


[نرن) ریت ررن6] <) پایه چرخه‌های کمینه نامیده 


می‌شود. در صورتی که مقدار زیر کمینه شود: 
(5) ,ط 

۱0-۶ ( 9 
1-21 


واضح است که (7)6(<1,6 پس یافتن یک پایه 
چرخه‌های کمینه متناظر با یافتن تنک‌ترین ماتریس ب) و 


از طرفی یافتن تنک‌ترین ماتریس ) می‌باشد. 
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پایه چرخه‌های شبه کمینه. پایه جرخه‌هایی که (1,0 
مقدار نزدیک به کمینه دارد. پایه چرخه‌های شبه کمینه 


نامیده می‌شود. 


پایه چرخه‌های بهینه. پایه چرخه‌هایی که منجر به تنکی 
ماتریس همسایگی پایه چرخه‌ها 6 می‌شوند. پایه 
چرخه‌های بهینه نامیده می‌شوند. بدیهی است چنین 
پایه‌هایی دارای کم‌ترین هم‌پوشانی هستند. 


پایه چرخه‌های شبه بهینه. اگر (7,)00 از مقدار مینیمم 
بیشتر باشد. چنین پایه‌ای را شبه بهینه می‌نامند. 
ضریب اشتراک ماتریس همسایگی سیکل عضو 
برای سطر نام برابر تعداد سطرهای [ام می‌باشد. به طوری 
که دو شرط زير برای سطر ام احراز شوند: 
(۱۵( ((٩),ط,‏ .... و 1+2 ,1+1 > [ 
شرط دوم به این معنی است که حداقل یک 5 پیدا 
شود که هر دو چرخه ز) و ز) شامل عضو 7 باشند و 
به بیان ماتریسی سطرهای ام و [ام ماتریس مجاورت 
چرخه عضو در ستون > درایه غیر صفر داشته باشند. 
به عنوان یک نتیجه ری ,ن) همواره برابر صفر است. نتیجه‌ی 
مستقیم تعریف اخیر به‌صورت رابطه‌ی زیر نشان داده 
می‌شود: 
1-(5) ,ط 
۱۹0( 4 2 ۱,۵2 (020) 
له (6)ن تشن خهه ی مسا روت هرس قی یا 
خودش می‌باشد. ضریب ۲ در جمله‌ی دوم به‌علت تقارن 
ماتریس "06 وارد شده است و حاصل جمع (6,)0 
نشان دهنده‌ی مجموع تعداد چرخه‌هایی است که با 
چرخه‌ی ام عضو مشترک دارند. پس 9 2 برابر تعداد 
2 
اعضایی غیر صفر و غیر قطری ماتریس )) خواهد بود. 
از آن‌جایی که (۳,)5 برای سازه‌ی معلوم مقدار معلومی 
می‌باشد. بنابراین تنک کردن ماتریس 06 تنها با مینیمم 
کردن مقدار رج) ی " !9" میسرخواهد بود یعنی چرخه‌ها 
طوری انتخانب ی از داشته باشند. 
روش‌های انتخاب چنین پایه چرخه‌هایی توسط کاوه 
[25-26-29-30-31] کتسن یافته است. 
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۳۱ 


با توجه به مطالب ذکر شده مناسب‌ترین نوع چرخه‌هایی 

که برای روش نیروها می‌توان انتخاب کرد چرخه‌های 

بهینه می‌باشند؛ ولی به دودلیل عمده‌ی زیر اکثر تحقیقات 

انجام شده بر روی چرخه‌های کمینه متمرکز شده است: 

۱. اکثر الگوریتم‌های انتخاب چرخه‌های بهینه پیچیده 
می‌باشند و از نظر زمان محاسباتی رایانه مطلوب 
تب باشتد. 

۲ باکر سازه‌های من شاب سیک هی که 


منجر به تشکیل چرخه‌های بهینه نیز می‌شود. 


انتخاب پایه چرخه‌ها 

پایه چرخه های یک گراف کاربرد زیادی در زمینه‌های 
گوناگون مهندسی دارد. میزان کارآیی هر یک از این 
کاربردها به پایه چرخه‌ی انتخاب شده بستگی دارد. یک 
پایه با چرخه‌های کوتاه‌تر زمان و حافظه‌ی مورد نیاز را 
رای ی از ها کاشش رات ده» شتن اشعاتایک 
پایه چرخه‌ی حداقل. مطلوب است و برای کاربردهای 
دیگری چرخه‌های با تقاطع حداقل مورد نیاز است؛ یعنی 
پایه چرخه‌های بهینه ترجیح داده می‌شود. پایه چرخه‌ی 
اقا هر انا بط اکتان ا ها مردد پروسی قران 
خرفت و ترس زیگرو] 13۵ اسلا گردیت: انگزریشهای 
عملی برای انتخاب پایه چرخه‌های نزدیک به حداقل 
توسط کاوه [11] و کسل وهم‌کاران [24] توسعه داده شد. 
روش‌های مشابهی برای تشکیل پایه چرخه حداقل یک 
گزافت فوشطیسار و هایکا |[ و2 اقا شیک فون 
مشابهی به‌وسیله‌ی کاوه [34] برای گراف‌های مستوی 
ساخته شد؛ با وجود این نمونه‌ی متقابلی هم برای 
گراف‌های فضایی داده شد. کولاسینسکو [35] یک نمونه‌ی 
روش ها سس متیر 
هورتون [36] یک الگوریتم از رتبه‌ی چند جمله‌ای برای 
یافتن پایه چرخه‌های حداقل گراف‌ها معرفی کرده است. 


تشکیل یک چرخه حدافل روی یک عضو. یک چرخه 
6 با طول حداقل روی یک عضو 70 که مولد نام دارد را 


می‌توان با استفاده از الگوریتم 7 کاوه به شرح زیر 
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تشکیل داد: 

دو 5۳1 از دو گره انتهای .0 و 7 عضو ۳0 تشکیل 
دهید و جریان را به برخورد دو 9187 با یکدیگر (عضو 
هه شامل 5107 نباشد) محدود کنید. جمم کوتاه‌ترین 
مسیر بین :2 و 16 :2 و 2 و عضو 1۳ یک چرخه حداقل 
ن) را روی 1۳ تشکیل می‌دهد. بااستفاده از ان 
الگوریتم‌ها طول‌ها را هم می‌توان محاسبه نمود. برای 
مثال. زنایک چرخه‌ی حداقل روی 7 می‌باشد؛ شکل(۱). 
7 ها با خط کلفت نمایش داده شده‌اند. به محض این 
که 1 پیدا شود مراحل تولید 1 قطع می‌گردد 


شکل(0:۱). 


به طور مشابه می‌توان یک چرخه‌ی حداقل» روی 
عضو 17 که از گره 1۷ عبور می‌کند. تولید کرد. یک 517 
که از :9 شروع می‌شود. تشکیل دهید و به محض این که 
به دو گروه انتهایی 170 رسید. جمع کوتاه‌تترین مسیرهای 
پیدا شده. بین ۰ 7 و 0 ۲ و ۲ با 1۳ جرخه‌ی مورد نظطر 
و تا مر دهد ای وود یرانق 
مشخص شلده است. 


ان فان مر هگا هر رها یک کنات :را 


می‌توان به عنوان بردار ستونی از ماتریس تلاقی جرخه - 
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استفاده از روش‌های نظریه گراف‌ها برای بهبود ... 


اعضاء مورد بررسی قرار داد. یک روش جبری, مانند 
حذف گوس را می‌توان برای وارسی یک چرخه نسبت به 
زیر- پایه انتخاب شده‌ی قبلی به‌کار برد. به هر حال. 
اگرچه این روش معمول است و ترتیب عدم استقلال 
الگوریتم‌های انتخاب چرخه را کاهش می‌دهد مانند 
خیلی از روش‌های جبری دیگر کاربرد آن نیاز به یک 
فضای ذخیره‌سازی قابل توجه دارد. 

اغلب روش‌های معمولی گراف تشوریکی از یک 
درخت گسترنده و استفاده می‌کنند و پایه چرخه‌های 
درختی آن را تشکیل می‌دهند. این روش خیلی ساده 
است؛ هر چند که معمولاً کاربرد آن چرخه‌های طولانی را 
نتیجه می‌دهد. این روش را می‌توان با اجازه دادن به 
انتخاب هر وتر به‌کار رفته در مجموعه شاخه‌های درخت 
انتخاب شده اصلاح کرد. تخفیف بیشتری در طول را 
می‌توان با تولید یک 511 از گره مرکزی گراف و به‌کار 
بردن این وترها به ترتیب صعودی فاصله‌ی آن‌ها از گره 
مرکزی به‌دست آورد. کاوه [11] 

روش سوم که گراف تثوریکی می‌باشد. شامل 
استفاده از چرخه‌های مجاز است. به فرآیند گسترشی زیر 
توجه کنید: 
642 یج راک رت را 

که > - *) یک چرخه بن) سیکل مجاز نامیده 
می‌شود. آ ال ببی‌لا؟) جب 6 داشته باشیم: 
)۱۸( 1+ (0 علاط < (6۳۳)رط 


به راحتی می‌توان اثبات کرد که شرایط مجاز بودن 
با ارضا یکی از شرایط زير تأمین می‌شود: 
ب ی که 0 یک تقاطع تهی است. 
۲ 2۲-5 (,ب,۵),ظ که ۲و 5 به ترتیب تعداد مولفه‌های 
وی 6 انتتن 
۳ 20 (بببه),ط که 0۳۳ و "6 پیوسته‌اند ۰( <6) 

در نسبت بالا 1+ 1,3 (5۸ که ,۷۲ و لا به 
ترتیب تعداد اعضاء و گره‌های :۸۵ می‌باشد. برای متنال» 
ترتیب انتخاب سیکل در شکل (۲): با شماره آن‌ها 


مشخص شده اقا 
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فرزاد شفیعی دیزجی- مصطفی حانزادی 


شکل۲ یک سیکل و سیکل‌های محصور آن 


روش متفاوتی توسط سیسلو وهوبی‌سکا پيشنهاد 
شده که: 
۵ (۵),<( )ره 
و برای مجازبودن مورداستفاده قرار می گیرد. به هر 
حال. یک پروسه تکمیلی برای 1 < « انجام می‌شود. برای 
مثال. وقتی که ون) به سیکل ۳) اضافه می‌شود اولین 
عدد بتی» کاوه [38] به اندازه‌ی ۳ افزایش پیدا می‌کند» 
بنابراین سیکل‌های ) و 2) قبل از افزايش انوی در این 
مرحله تولید می‌شود. 
در این مقاله برای وارسی استقلال چرخه‌ها از 


روش اول وسوم استفاده شده است. 


پایه چرخه‌هایی با شرایط موزونی بهینه 

به منظور بهینه کردن شرایط موزونی ماتریس‌های نرمی» 
پایه‌های استاتیکی مربوط به پایه چرخه‌هایی که خواص 
ویژه‌ای دارند باید انتخاب شود. یک پایه چرخه اگر 
خواص زیر را داشته باشد به‌عنوان یک پایه چرخه با 
شرایط موزونی بهینه تعریف می‌شود: 
(الف): پایه چرخه بهینه باشد, یعنی عناصر غیر صفر 
مربوط به ماتریس مجاورت چرخه‌ها حداقل باشد که به 
ی تاک تفت جر ور 
(ب): اعضاء مشترک چرخه‌ها با بیشترین وزن باشند؛ 
یعنی عنصرهای غیرقطری مربوط به ماتریس نرمی 
کوچکترین مقدار بزرگی ممکن را داشته باشند. 

یک گراف وزن‌دار ممکن است بیشتر از یک پایه 
چرخه‌ی بهینه داشته باشد. ارضا شرط(ب) شرایط 


موزونی را بهینه می‌کند. به هر حال, اگر یک چنین پایه 
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شرط‌های (الف)و(ب) پیدا شود. به‌عبارت دیگر پایه‌ای 
را باید انتخاب نمود که به‌طور جزئی هر دو شرط بالا را 
ارضا کند. از آن‌جایی که الگوریتمی برای تشکیل پایه 
چرخه بهینه موجود نیست. تحقیق برای پایه چرخه‌هایی 
با شرایط موزونی نزدیک به بهینه انجام می‌شود. مسأله‌ی 
انتخاب پایه چرخه‌ی بهینه را می‌توان با شکل ریاضی زیر 


توضیح داد: 
۲ 1-(6) رد 
1۹ 10۲ ۸:8 
1-1[ 
0 1-(5) ,۲ 
۲۰( (رب6] ۷۷۵ ۱۵ 
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که 5 انقباض یافته‌ی 5 است و ,6 [ا< 0 مي‌باشد. 
از 
می‌توان یک که میا له ایک عسا لش بهینه‌سازی جند 
موضوعه است و الگوریتم‌هایی که در ادامه می‌آبد» به 


گونه‌ای تعریف شده‌اند که هر دو تابع موضوع را به طور 
سر مهیتشانرعیا خن کاود [37] 


خطاصا. تحلیل خطا در یک محاسبه‌ی عددی. 
به ندرت دقیق هستند. چرا که آن‌ها در اکثر مواقع بر 
پایه‌ی مشاهدات و برآوردها بنا شده‌اند و از طرفی, حود 
روش‌های عددی. موجب انواع مختلفی از خطاها 
زير به ذکر اهمیت آن‌ها می‌پردازيم: 

فرض کنید بخضواهيم ریشه‌ی معادله‌ی 
0- 0.00008 + :0.4002 + 7« راباع رقم اعشار 
حساب کنیم. جواب به دست ا فتاه برابر ۳ 

جواب حقیقی به‌دست آمده با دقت ۸ رقم اعشان 
براپر ۰/۰۰۰۲- می‌باشد؛ بنابراین حساب با اعداد دارای ۶ 
پدیده‌ی مزبور استفاده از حساب با اعداد اعشاری ۶ 


رقمی می‌باشد. اما نباید تصور کنیم که اعداد اعشاری ۸ 
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۳ 


رقمی تمام مسائل را حل می‌کنند. 
حساسیت (بد وضعی) و ناپایداری. برخی محاسبات 
به خطای گرد کردن بسیار حساس می‌باشند؛ در حالی که 
برخی دیگر این گونه نیستند. در پاره‌ای موارد می‌توان 
این حساسیت را با تغییر روش محاسبه از بین برد لیکن 
در بسیاری از موارد مسائل به طور ذاتی به خطای گرد 
کردن حساس می‌باشد. بنابراین باید حساسیت به روش 
و حساسیت ذاتی را از یکدیگر متمایز ساخت.کلمه‌ی 
حساسیت به طور کلی به وابستگی پیوسته‌ی یک مسأله 
نسبت به داده‌های ورودی و يا روش آن اطلاق می‌گردد. 
چنانچه خطا در یک مسأله با نرخ چند برابر مرتب 
افزایش یابد. روش را از لحاظ عددی, ناپایدار می‌نامند. 
مشال: معادله‌ی مرتبه چهار 2-0 (6-2 را از نطر 
بدحالتی بررسی کنید. 
حل: این معادله دارای چهار ریشه‌ی یکسان 22 می‌باشد. 
ال ای ماوت وان تک عوو سای رک رما تنیز 10۳ 
کسر و رفتار آن را مورد مطالعه قرار می‌دهیم: 
(0۲۱ . 1599999999-0+-160- ۵+ مد << 0-2۳-10 
ریشه‌های این معادله باید رابطه‌ی *10+< (2-) 
را ارضا نمایند یعنی *4/210-(6-2. بنابراین 
ریشه‌های معادله عبارتند از: 
2-11 2 5 بف22+0.01 و5 ,21.99 بت ی01. 2 رتز 
همان‌گونه که ملاحظه می‌گردد. ریشه‌های معادله 
دافهشله میت به-داده‌ها قوال از قبیل یک ضریب 
بسیار حساس می‌باشد. این نوع رفتار تابع که مستقل از 
روش حل آن بوده و جزو خحصوصیات ذاتی آن می‌باشد» 
منجر به بدحالتی سیستم می‌گردد؛ چنانچه در این‌جا 
ریشه‌های معادله با اندکی تغییر در ضریب ثابت آن, دچار 
تغییرات بسیار بزرگی گردیده‌اند. به‌طورکلی, یک ماتریس 
اک تطداقل یکی از فشته های. زیر رادشه باشعکه 
بد وضع نامیده می‌شود: 
تغییر کوچکی در ضرایب (عناصر ماتریس) باعث 
تغییرات بزرگی در جواب دستگاه شود. 
۲ عناصر قطری ماتریس, نسبت به ساير عناصر آن 
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کوچکتر باشند. 
۳ مقدار | ۸۵۱۰/۸ دقیقاً یک نشود. 
6 مقداز "1" ه] دقیقا با #برایر باشد: 
۵ مقدار [*۵۵] دقیقاً با ماتریس یکه آ برابر نباشد. 
1 مقدار [*[(8)] ۸] از ماتریس یکه آ بیشتر از 


مقدار ۸۵ از [ تفاوت داشته باشد. 

ی مه ی رش ان کم اش 
محورگیری بر افزايش دقت جواب تأثیر دار اما در مورد 
مسائلی که بد وضعی شدیدتری دارند فقط با دقت 
اشتا سا سای مر ترا ادفت میات( آفرایش فاد 


اعداد شرایط موزونی. به منظور اندازه‌گیری شرایط 
موزونی یک ماتریس, اعداد گوناگونی تعریف شده و به 
کار می‌رود. سه عدد شرایط موزونی که استفاده‌ی آن‌ها 
آسان است در ذیل شرح داده می‌شود؛ کاوه [38]: 
الف) نسبت مقادیر ویژه‌ی حداکثر و حداقل: مقادیر ویژه 
و بردارهای ویژه به شرایط موزونی ماتریس‌ها مرتبط 
است. نسبت مقادیر ویژه‌ی حداکثر و حداقل یک ماتریس 
ایس(۱/۱ سس را می‌توان به عنوان یک عدد شرایط 
موزونی به کار برد. به راحتی می‌توان نشان داد که 
لگاریتم در مبنای ۱۰ این عدد. شرایط موزونی به طور 
قابل توجهی با تعداد حداکثر ارقامی که در معکوس کردن 
یا حل معادلات هم‌زمان از دست می‌رود. متناسب است. 
بنابراین تعداد ارقام خوب ی در حل معادله برابر است با: 
)۳۲( 2۲-۲ ( مرس۱2/ ۱ )۲-108 2 8 
در این نسسبت 2( ۱/۱ 1080۸ باظر2 
عددی است که از یک ماشین به ماشین دیگر تغیسر 
می‌کند. برای مثال در 191/360 تقریباً ۸ رقم برای دقت 
ساده و ۱٩‏ رقم برای دقت مضاعف می‌باشد. باید ذکر 
شود که برآورد بالا محافظه کارانه است و تجربه نشان 
داده است که بآ یک رقم در جهت اطمینان است. 
ب) دترمینان ماتریس نرمالیزه شده‌ی سطری: یک مقدار 
کارا و ساده که شرایط موزونی یک مجموعه معادلات را 
ارزیابی می‌کند. دترمینان ماتریس نرمالیزه شده‌ی سطری 
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ماتریس ») , :ت) تقسیم می‌شود به: 
2 


۲۳( [ و بر در + ۳ 
بزرگی دترمینان ماتریس نرمالیزه شده سطری که با 
۷ مشخص می‌شود. معیار خوبی برای سنجش شرایط 
موزونی ماتریس 6 می‌باشد. واضح است که بزرگی ایسن 
دترمینان در فاصله‌ی زیر قرار دارد: 
۲( 1 ۲( > 0 
از آن‌جایی‌که 6 لزوماً معین مثبت است. در ماتریس 
با شرایط موزونی ایده آل ۳3-1 می‌باشد. این مورد 
زمانی اتفاق می‌افتد که ماتریس متعامد با قطری باشد؛ 
کاوه[38] . 
چ‌( نسبت دترمینان‌ها: از آن‌جایی که بهترین شرایط 
موزونی برای معکوس کردن در یک ماتریس قطری بودن 
آن می‌باشد. پارامتر زیر را می‌توان به عنوان یک معنی 
کاربردی برای اندازه‌گیری شرایط موزونی ماتریس 6) 
قبول کرد: 
(۲۵( [, 06۱01-16110 جع 
که رت)؛ ماتریس قطری مرکب از درایه‌های قطر 
ماتریس ت) می‌باشد واع0 ؛ معنی دترمینان را می‌دهد. 
مقادیر 6 نزدیک صفر برای یک ماتریس با شرایط 
موزونی ایده‌آل می‌باشد. بنابراین. عدد شرایط موزونی 
1 به شرح زیر استفاده می‌شود: 
۲0( [, 16116 /[06۱]01 - 2۳۲۲ 
برای یک ماتریس با شرایط موزونی ایده آل ۳۳۳۲ 
نزدیک واحد است؛ کاوه [38]. 


گراف وزن دار و عضو محاز. سختی (يا نرمی) نسبی 
اعضاء یک سازه را می‌توان به عنوان اعداد صحیح مثبت 
مربوط به اعضاء گراف مدل یک سازه مورد بررسی قرار 
داد که یک گراف وزن دار را نتیجه می‌دهد؛ کاوه [38]. 
5 را مدل یک سازه‌ی قابی در نظر بگیرید و زک به 
اصلی انتخاب له برای سازه تعریف می‌شود. وزنی را 
برای :1۳ می‌توان با استفاده از عنصرهای قطری تیک از نم 
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تعریف کرد: 
(0 + 0 ۳ 2601 "1 کر < ( ,)۲۱۷ 
[ )+ یم + یم 2 << (ر) ۱۷ 
۳۳( 


8۸ 1۳1 
0, < , 0 < 


0 69( 
۳ 


1 ۳ 
در تعریف دیگری می‌توان مجموع ریشه‌ی دوم 
عنصرهای قطری مب را به کار برد.توابع دیگری نیز برای 
تخایکی تن نس ایام ۵ من وان تعرفت کر 
تعریف: فرض کنید وزن اعضاء ۳ ر ۳02 ر ...۰ , (1 به 
ترتیب با (,۳) ۷۷( :)۷۷ .۲۷۷)(:۰ تعری_ ف 
می‌شود. عضو ۲ یک عضو 20۳551016 -۲) نامیده 
می‌شود اگر: 
۳۹۱( ۱ 


1/6( 
[1 


۷۳ 
0 
که 0 یک عدد صحیح است و می‌تواند ۲ ۲ و ... 
فرض شود که ما 0 را برابر ۲ در نظر می‌گيريم. به هر 
حال. یک مطالعه‌ی کامل برای مقادیر دیگر 4 لازم است. 
اگر یک عضو مجاز نباشد. غیر مجاز یا 5- مجاز 


و20 -5) نامیده می‌شود. 


تعاریفی از تتوری گراف‌ها 
یک گراف شامل یک مجموعه (68 از اجزایی موسوم 
به گره (نقطه‌ها یا رأس‌ها) و یک مجموعه (11)5 از 
اجزایی موسوم به اعضا (پاره خط‌ها یا منحنی‌ها) است که 
هر عضو با یک جفت گره که انتهای آن نامیده می‌شوده 
مربوط است؛ کاوه [38]. 

اگر دو يا چند عضو متصل به گره‌های یکسان 
باشند به آن‌ها اعضا مضرب می‌گویيم و یک عضوی که 
دو انتهای آن به یک گره متصل باشد. حلقه نامیده 
مین ود 

گرافی که حلقه و اعضاء مضرب ندارد گراف 
ساده نامیده می‌شود. اگر (1168 و (1۷1)5 قابل شمارش 
باشند. گراف مربوطه محدود نام می‌گیرد. گراف‌های 
به‌کار رفته به‌عنوان مدل ریاضی سیستم در سازه‌ها گراف 
ساده محدود می‌باشند. دو گره از یک گراف را همسایه 
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می‌نامند. هرگاه این گره‌ها در دو انتهای یک عضو واقع 
شده باشند. یک عضو متلاقی با یک گره خواهد بود. اگر 
آن گره. یکی از گره‌های انتهای عضو باشد. یک زیرگراف 
ز5 از ٩‏ خود گرافی است که (1168> (:1166 ,(1۷65- ( 6 
است و گره‌های انتهایی هر عضو از :5 همان گره‌های 
انتهایی در گراف ٩‏ می‌باشد. یک گام ۳ از 5 یک رشته‌ی 
محدود و وک ۳ نت کته 
جمله‌های آن متناوباً گره‌های :۲ و اعضا بط از ٩‏ برای 
0 > 161 است و ۲ و 1 دو انتهای :10 می‌باشد. 

مسیر با عضو تکراری گامی است که هیچ عضوی 
از 5 بیشتر از یک بار در آن ظاهر نشود. طول مسیر :۳ با 
((1.)۳ نشان داده می‌شود و برابر تعداد اعضای آن است. 
7 کوتاه‌ترین مسیر بین دو گره :10 و ما است. اگر برای 
هر مسیر ۳ دیگری بین این گره‌ها (1.)۳(>,۳ باشد. 
فاصله بین دو گره از یک گراف تعداد اعضا کوتاه‌ترین 
مسیر بین اين دو گره است. اگر یک مسیری بین آن‌ها 
موجود باشد. به دو گره ز و زامتصل گفته می‌شود. 
گراف ٩‏ اگر همه جفت گره‌های آن به هم متصل باشد 
گراف یکپارچه و در غیر این صورت گراف چند پارچه 
نامیده می‌شود. 

تعداد پارچه‌های گراف ٩‏ حداکثر زیرگراف‌های 
بکتارشته آن ا دک پم کف نیع اند سیخ 
( 10 و133 ود وی وبثلا ر0) که بظع 1 و 821 می‌باشد؛ 
یعنی چرخه یک مسیر بسته است. به طور مشابه یک 
چرخه با گره تکراری بسته (چرخه لولادار) و یک گام 
بسته را می‌توان تعریف کرد. درخت ]از 5 یک 
زیرگراف متصل از 5 است که شامل چرخه نباشد. یک 
مجموعه از درختان 5 یک کل .وا تشکیا: شي‌دهنت: اگر 
یک درخت. شامل تمامی گره‌های 5 باشد. درخت 
گسترنده نامیده می‌شود. برای سادگی به آن به‌عنوان یک 
درخت اشاره می‌شود. درخت کوتاه‌ترین مسیر 91 از 
یک گره مشخص 10 ريشه می‌گیرد و درحتی است که 
فاصله‌ی بین هر دو گره :7 و 70 کمینه است. دو فضای 


برداری به گراف 5 ارتباط داده می‌شود. این فضاها 
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عبارتند از: فضای چرخه و فضای مجموعه برشی. 
با این که فضاهای برداری را می‌توان روی یک 
میدان دلخواه تشکیل داد. برای سادگی. میدان از مدول 


صحیح ۲ مورد بررسی قرار می‌گیرد که در آن 1+1-0. 


چرخه‌ای را بررسی می‌کنیم. مجموعه‌ی تهی نیز به‌عنوان 
یک عضو در نظر گرفته می‌شود. می‌توان نشان داد که 
بردار مجموعه چرخه‌ای است. بنابراین بردارهای مجموعه 
چرخه‌ای» یک فضای برداری می‌باشد که بعد آن برابر 


(۳۰( (5),ط + (۲۲)6 - (۳166 - (5) ,ط > (۲۵01>)6 


پایه های جرخه‌ای اساسسی. برای فضای مجموعه 
چرخه‌ای می‌توان یک پایه تشکیل داد. به درخت 1و 
همراه درخت "۲ توجه کنید. زیرگراف 5 که مرکب از 
همراه درخت‌ها می‌باشد و هر عضو "1 (وتر) یک 
مجموعه چرخه‌ای که به‌عنوان مجموعه چرخه‌ای اساسی 
شناخته می‌شود را تشکیل می‌دهد. مجموعه‌ای از 
مجموعه‌های چرخه‌ای که با اضافه کردن شاخه‌های "۲ 
به 1 تهیه می‌شوند. با هم پایه‌ای را برای فضای مجموعه 
چرخه‌ای 5 تشکیل می‌دهد که به‌عنوان پایه چرخه‌ای 
ات رم شوه 


ماتریس‌های پایه‌های چرخه‌ای. ماتریس تلاقی عضو 
مجموعه چرخه‌ای» ) برای یک گراف» 5 یک سطر به 
ازای هر مجموعه چرخه‌ای ٩‏ و یک ستون به ازاء هر 
عضو دارد. درایه رن) از ن) برابر ۱ است. اگر مجموعه 
چرخه‌ای ,) شامل عضو 50 باشد و در غیر این صورت 
برابر صفر است. سطرهای مستقل :) برای یک پایه 
برشی با ن) نشان داده می‌شود و ماتریسی با عنوان 
ماتریس تلاقی پایه چرخه‌ای تشکیل می‌دهد که یک 
ماتریس 5(*1۷), است. ((۲,)5 رتبه گراف را نشان 


می‌دهد). ماتریس همسایگی مجموعه چرخه‌ای 
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2-00 ظ یک ماتریس (8(*0,)5), است اگر با و 
رن) یک عضو مشترک داشته باشند. هر درایه رنه آن برابر 


۱ است؛ و درغیر این‌صورت. برابر صفر است؛کاوه [38]. 


الگوریتم‌هایی برای تشکیل پایه چرخه‌های نزدیک 
به بهینه 

در این بخش دو الگوریتم برای انتخاب پایه چرخه‌هایی 
با شرایط موزونی نزدیک به بهینه یک گراف وزن‌دار 
شرح داده می‌شود. در هر چرخه‌ی انتخاب شده. سه یا 
شش (8.۳.59) مربوط به 5 که به ترتیب از بک قاب 
مستوی پا فضایی می‌باشد تن نشکا می‌گردد؛ کاوه 
[37-38]. 


الگوریتم اول. مرحله‌ی ۱: از یک عضو گراف یک 
چرخه‌ای با طول مینیمم و وزن ماکزیمم ایجاد می‌کنیم. 
ریمعت ال تزانی ق تون )افیا رازن 
انجام می‌دهیم؛ یعنی از تمامی اعضای گراف مورد نظر 
چرخه‌ای می‌اندازيم که در حین مینیمم بودن طول چرخه 
ماکزیمم وزن را نیز داشته باشد. 

(از دو انتهای عضو 91 باوزن ماکزیمم 
به‌صورت هم‌زمان تولید می‌کنیم و به محض این که دو 
۲ با ماکزیمم وزن همدیگر را قطع می‌کنند با یک 
مسیر (۳۵) به دو انتهای عضو مورد نظر برمی‌گردیم. 
عضو مورد نظر و دو مسیرایجاد شده تشکیل چرخحه 
مینیمم را می‌دهند؛ بدین صورت. یک چرخه با طول 
مینیمم و وزن ماکزیمم تولید می‌کنیم.) 
طرز تشکیل 911 ماگزیمم به‌صورت زیر است: 
الف) تمام اعضایی که با 0 تلاقی دارند را در نظر بگیرید 
تیا رای ای سگوسان یبای نها 
ب) تمام اعضایی که با گره‌های مشخص شده با عدد ۱ 
تلاقی دارند را بيابید و آن‌ها را به طور صعودی بر حسب 
وزنشان مرتب کنید. 
ج) اعضاء درختی را از اعضاء مرتب شده‌ی بالا انتضاب 


کنید و عدد ۲ را برای انتهای دیگرشان اختصاص دهید. 
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مرحله‌ی (ج) را تا زمانی ادامه می‌دهیم که دو 5117 
همدیگر را قطع کنند. 
مرحله ۳: بعد از اين‌که مرحله‌ی ۲ را به اتمام رساندیم و 
از تمام اعضاء گراف یک چرخه با طول مینیمم و وزن 
ماکزیمم تولید کردیم وزن هر چرخه را به دست 
می‌آوریم. (وزن یک چرخه به‌صورت مجموع وزن 
اعضای آن به دست می‌آید.) 
مرحله‌ی : چرخه‌های به دست آمده را به‌صورت 
صعودی وزن در یک مجموعه مرتب (0۳067) می‌کنیم. در 
واقع جرخه‌های به‌دست آمده را 50۲0) می‌کنیم. 
مرحله‌ی ۵: از الگوریتم گریدی (ط۸۱۵0۲ 60ع0۳) 
استفاده می‌کنيم و چرخه‌های مورد نیاز را تا وفتی که 
(۲:66 تا چرخه‌ی مستقل به دست آید, انتخاب می‌کنیم. 
(از الگوریتم گریدی ابتدا یک چرخه با ماگزیمم وزن از 
مجموعه مرتب (0۳06۲) شده‌ی مرحله‌ی ۶ انتخاب 
می‌کنیم و سپس چرخه‌ی بعدی با وزن ماگزیمم از 
مجموعه‌ی باقی‌مانده انتضاب می‌کنیم که مستقل از 
چرخه‌ی به‌دست آمده‌ی اولی است و مرحله‌ی پنجم را 
شانوهانی که (یط تا چرخهی تفکیل رد ادا 
می‌دهيم. بدین ترتیب یک پایه‌ی چرخه با شرایط 
موزونی نزدیک به بهینه تشکیل می‌گردد.) 


الگوریتم دوم. مرحله‌ی ۱: مدل گرافی سازه‌ی مورد 
نظر راپیدا کرده. سپس از یک عضو گراف مورد نظ 
چرخه‌ای با وزن ماگزیمم وطول مینیمم تولید می‌کنیم. 

ی ۱ هرارق را ای ای اعشای اف 
انجام می‌دهیم؛ یعنی از تمام اعضای گراف مورد نظر 
چرخه‌ای می‌اندازيم که در عين مینیمم بودن طول چرخه. 
ماگزیمم وزن را نیز داشته باشد. البته منظور ما از جمله‌ی 
بالا این است که بین مینیمم شدن طول چرخه و ماگزیمم 
بودن وزن چرخه سازشی ایجاد کنیم؛ زیرا نمی شود هم 
طول چرخه را مینیمم و هم وزن چرخه را ماگزیمم کرد. 
همان طور که قبلا هم اشاره شد. پیدا کردن الگوریتمی که 
پرصفرترین وخوش حالت ترین ماتریس را برای ما تولید 
کل ورد ردو نت کل ان بخ مس نا روت 
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سازی چند منظوره است و باید بین این دوخاصیت 
سازشی ایجاد کنیم؛ یعنی اگر در مساله‌ای پرصفری برای 
ما حائز اهمیت باشد از الگوریتم‌هایی استفاده کنیم که 
ماتریس ما را پر صفر می‌کند و برعکس این مسأله هم 
صادق است. (از دو انتهای عضو 9۳1۷ بهبود یافته با 
وزن ماکزیمم به‌صورت هم‌زمان تولید می‌کنیم و به 
محض این که دو 91*11۷ بهبود یافته با ماکزیمم وزن 
همدیگر را قطع می‌کنند. با یک مسیر (۴۵1۳) به دو انتهای 
عضو مورد نظر برمی گردیم. عضو مورد نظر و دو مسیر 
ایجاد شده. تشکیل چرخه‌ی مینیمم را می‌دهند؛ بدین 
صورت. یک چرخه با طول مینیمم و وزن ماکزیمم تولید 
طرز تشکیل ۹۳۲1۷ بهبود یافته به‌صورت زیر است: 
گرهی از آن‌جا که 98171 را قرار است شروع 
کنیم, مشخص کرده و به آن ٩0۳‏ می‌گذاريم و سپس 
عضوهایی را که با "۲0 تلاقی دارند مشخص کرده و 
سپس متوسط وزن عضوهای متقاطع با "0" را محاسبه 
عضوهاست (عضوهای متقاطع با 0۳؟) از شاخه‌های 
درحت 7 شدن محروم می‌کنيم. در ادامه تمام 
اعضایی را که با ۵ تلاقی دارند در نظر بگیرید و عدد ۱ را 
برای انتهای دیگرشان اختصاص دهید. سپس تمام 
اعضایی را که با گره‌های مشخص شده با عدد ۱ تلاقی 
دارند بیابید و آن‌ها را به طور صعودی بر حسب وزنشان 
مرتب کنید. در ادامه اعضاء درختی را از اعضای مرتب 
شده بالا انتخاب کرده و عدد ۲ را برای انتهای دیگرشان 
احتصاص دهید. سپس برای گره‌هایی که با عدد ۱ 
مشخحص شده است نیز متوسط وزن اعضای درختی بالا 
را به دست می‌آوریم و اعضایی که وزنشان از وزن 
مترشط نهدست آمده کم است: زا از 5911 شبن 
حذف می‌کنیم. سپس این عملیات را برای گره ۲ و سایر 
گره‌های گراف انجام می‌دهیم تا این که 1۷] 9۳ ما تشکیل 
گرگ 
مرحله‌ی ۳: مرحله‌ی دوم را برای تمام اعضای گراف 
مینیمم و وزن ماگزیمم تولید کردیم. 
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استفاده از روش‌های نظریه گراف‌ها برای بهبود ... 


سپس وزن هر چرخه را به دست می‌آوريم. (وزن 
یک چرخه به‌صورت مجموع وزن اعضاء آن به دست 
می‌اید.) 
مرحله‌ی ۶: چرخه‌های مخت امسژه را به‌صورت 
صعودی وزن در یک مجموعه مرتب می‌کنیم. 
مرحله‌ی ۵: از الگوریتم گریدی (۸۱۵0۲۳ 66ع0۳) 
استفاده می‌کنيم و چرخه‌های مورد نیاز را تا وفتی که 
(5),ظ تا چرخه مستقل به دست آیده انتخاب می‌کنیم. (از 
الگوریتم گریدی ابتدا یک چرخه با ماگزیمم وزن از 
مجموعه مرتب (0۲06۲) شده‌ی مرحله‌ی ۶ انتخاب 
می‌کنیم و سپس چرخه‌ی بعدی با وزن ماگزیمم از 
مجموعه‌ی باقی‌مانده انتخاب می‌کنیم که مستقل از چرخه 
به دست آمده‌ی اولی است ومرحله‌ی پنجم را تا زمانی 
که (9:)5 تا چرخه تشکیل شود ادامه می‌دهيم. بدین 
ترتیب» یک پایه چرخه با شرایط موزونی نزدیک به بهینه 


تشکیل مین کر کی 


نمونه مثال ها 
به منظور بررسی کارایی الگوریتم‌های ارائه شده برای 
الگوریتم‌های اول ودوم ارائه شده با زبان برنامه نویسی 
جاوا (18۵7۷۸۵) برنامه به‌صورت ویژوال نوشته شده است 
وخوش حالتی این دو الگوریتم با یکدیگر مقایسه شده 
است. شکل ویژوال برنامه به‌عنوان نمونه در شکل (۲) 
نمایش داده شده است. 
مشال ۱: یک قاب ۳ طبقه‌ی ۶ دهانه در شکل (1) 
ملاحظه شده است. سه حالت با استفاده از دو نوع 
خاصیت عضو مطالعه می‌گردد: 
هلا 521۷105 *۳ 0,00000171<,] *۳ 0.00106عب۸۵ 
گصل 7 10 2.17 هه 0:00019610< 1 5 0.00970عي۸ 
برای تمامی عضوهای سازه‌ی طول. یعنی 1,237 

قرار داده شده اشت: 
مثال ۲: یک قاب ۳ طبقه‌ی ۳ دهانه در شکل (۷) ملاحظه 
شده است. سه حالت با استفاده از دو نوع خاصیت عضو 
مطالعه یی که 

کل 10 52.17 *۳ 20:00000171< ,1 ۳ 0,00106-ر۸۵ 


کل 7 10 «2.1عظ هط 0,00019610< 1 0 0.00970< ر ۸ 


سال بیست و دوم شماره یک» ۱۳۸۹ 


فرزاد شفیعی دیزجی - مصطفی حانزادی ۳۹ 


با 


وا ام فمیضا 


مسا اس دوف سول حصه ها | «قمه: فت | ۵۵ سل [ آ] 
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شکل ۳ شمای برنامه‌ی به ویژوال در آمده 
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69 )0( 010( 
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شکل 1 ماتریس‌های هم الگوی ماتریس نرمی سازه برای شکل (4) 
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استفاده از روش‌های نظریه کراف‌ها برای بهبود ... 


جدول ۱ اعداد شرایط موزونی با الگوریتم اول 


۳۳ 
7.134199۳-0 
0 
0 


۷ 
6.381245۳-9 
0 
0 


بآ ۷۳ 1 
34 )2( 
145( )5( 
۹+6 0( 


جدول ۲ اعداد شرایط موزونی بهبود يافته با الگوریتم دوم 


۳۳ 


7.134199۳-0 
3.613358۳-1 
3./737124۳-30 


(12 


۳ 


6.381245۳-9 
2.06309 ۳-2 
8.307591۳-1 


(ظ) 


(6) 


م۳1 ۷۲ 
4 3 )0( 
4-18 [3) 
3.45 63 

۳9( 


0) 
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1 عا 2 
س زد عد 2 
۴ زد 
ت‌ له ت 
3 اد عد عد ت 
7 اد ود ت 
4 اد 
اد لد 2 ی 
. زد س 


شکل ٩‏ ماتریس‌های هم الگوی ماتریس نرمی سازه برای شکل (۸) 
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فرزاد شفیعی دیزجی- مصطفی حانزادی 


جدول ۳ اعداد شرایط موزونی با الگوریتم اول 


۱ 


۱92۵ ۳1 7 (1 26) 
(0) 2321214 5.7325461۳-1 2.214527۳5-4 33 
0) ۹-48 0 0 33 
9( (۹-6 691879۳5-0 3.524312۳8-1 33 

جدول؛ اعداد شرایط موزونی بهبود يافته با الگوریتم دوم 
(26)0 ۳ 7 1 ۳۷۳ 

33 22145275-4 5.7325461۳-1 321214 1 
45 45725415-6 32452487-4 58 0 
42 1.6848295۳-1 19220135-0 9 32( 3 
۳ ۰ اکسع 4 ىَ 

برای تمامی عضوهای سازه‌ی طول» یعتسی 1-۳[ نتبحه گیری 


راز داد بای 

برای مثال اول» چرخه‌های به‌دست آمده از 
الگوریتم اول در شکل (۵) و همچنین چرخه‌های به‌دست 
آمده از الگوریتم دوم در شکل (4) و ماتریس‌های هم 
الگوی ماتریس نرمی سازه با الگوریتم دوم در شکل (0) 
نشان داده شده‌اند. هم‌چنین در جداول (۱) و (۲) نتایج 
حاصل از الگوریتم‌ها به منظور سنجش خوش وضعی 
فاتر ین هار مسازهات دادم شنوانه: 

برای مشال دوم نیز مانند مشال اول. چرخه‌های 
به‌دست آمده از الگوریتم اول در شکل(۷) وهم‌چنین 
چرخه‌های به‌دست آمده از الگوریتم دوم ارائه شده در 
شکل(۸) و ماتریس‌های هم الگوی ماتریس نرمی سازه 
با الگوریتم دوم در شکل )٩(‏ نشان داده شده‌اند. هم‌چنین 
در جداول (۳) و (6) نتایج حاصل از الگوریتم‌ها به 
منظور سنجش خوش وضعی ماتریس‌های سازه‌ای داده 


شده‌اند. 


سال بیست و دوم شماره یک. ۱۳۸۹ 


برای هر سازه رشته‌ای می‌توان یک مدل کر اقنتي انتخاب 
کرد و مسأله‌ی تحلیل بهینه سازه‌ای اک میا 
بهینه‌سازی ترکیباتی تبدیل کرد. هم الگویی بین 
ماترس‌های سازه‌ای و ماتریس‌هایی که با مدل گرافی 
سازه هم‌پست شده‌اند اساس اصلی تخییر مسائل تحلیل 
بهینه سازه‌ای به مسائل بهینه سازی توپولوژیکی است. 
درمثال‌های ارائه شده نشان داده شده است که می‌توان با 
ماتریس‌های سازه‌ای را بهبود بخشید. در مثال‌ها دیده شد 
که خوش حالت شدن ماتریس نرمی, باعث به وجود 
آمدن غیرصفرهایی در ماتریس هم الگو با ماتریس نرمی 
و در نهایت خود ماتریس نرمی می‌گردد و اين جریمهای 
ماتریس نرمی و هم‌زمان پرصفر کردن ماتریس نرمی. یک 
مسأله‌ی بهینه سازی چند منظوره است و باید سازشی بین 


این دو خاصیت ایجاد شود. 
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۶:۲ استفاده از روش‌های نظریه گراف‌ها برای بهبود ... 
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0۶۵ عا معصنعدا ععتتت‌نتتای تصعصعاه . ماتص ۵۴ وتولقمه »ءتاقصعافوه معط ما ممتاتاطتتاجی ۳۶۸ رظ ب0۳]" 
1979(۰) ,.01۷76۲] 56و رات6ککن0 100۵0۲۵1 ررصقصصهی مق 2۵00۳ 
بل ر"عمطاعص م۵۲ عطا ۵ ومت0یميم امصمتامانامم ‏ ۵ ری رگناتتعنط] 20 1۷۲ ر1۵2۷0۵ ربا و1۵۳6 
1469-1495)1982(۰ ۲0۰ ,18 .۲۱۵08 ۱۷۲6۰ .تعصصت ۲ 
و.۷۲6۸ 5۰( .و۸۱ .[ 51۸0۷ ,"ع)ووولاناط معط ۶۵۲ وتعودا ععتممو ج هصتانام‌ومن .۲ ,۱ رطلهع۳۱ 220 ربط.ژ راتعطاا6 
446-459)1987(۰ .00 ر8 
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فرزاد شفیعی دیزجی- مصطفی حانزادی 


۰ ,35 یکلتاقع ۱۷۱۵ عط)وزمصص رصمتاهمتصتامن امست‌بوی صا ممتاه)نامرصوی فتقفقط ات 592786 ریق رحعطا0ظ 
501-519)1989(۰ 

و7 و ۷/۵۰ 56 .و2 .[ 54 روتها۵(۵ ۲۱ صماطامتج مهو ان عظ۲ ربظ محعطا0ظ 220 ر.۲ .۲ رصقفصامن 
۰ (527-5337)1986 2۳۰ 

و8 ,۷/6 156۰ .۸9 . 54 ,هامرهم 11 جماطامتم معهوو ((تاظ م۳۲ ربظ مطعطا۵ 220 .۲ .] مرصقصامن 
(544-561)1987 .0۳ 

ح1ح1 ,00 ر22 رتاک تامصصمت ر0مطاعص م۵۲ 0تجصهای عمط عناعیع۷ 0مطاعص م۲۵۲6 ۰۱۳۱6۵۲۵۲60 ر.ل .5 رعالظ۳۵0 
164)1986(۰ 

-129 .09 ,24 .[ ۸744 ر"00اع م1۵۲6 0عاجععاص عطا گ۵۶ صمتاعاصم تمصمتاعتنه۷ م۲ رل .5 رکالط/۳2 
13761986 

.0 .ل ,"۱961000 10۲6۵ 0عاهتععاصا رها عمستم‌نتای ۵ فتوررآمهه ۲۳۲۵0۵۵ ر.و رتتتعل۷2 280 .لا .5 کانقط۳20 
(93-109)1982 ۲۳0۰ ر83 و.1۷10 

۰ 0130 طا وععط۸072 )۱6 رت طججم ع ۵۲ ععلبن 0۲ وعووط امصصتصن۱ رب۷۲ رملورک ۵ص .۲ رقیازطانا۲۲ 
(283-293)1975 0۰ ,۲۲665 ۸6۵06۲۲6 

۶ وزوت‌امصم نازدانه عمط >م۶ وععودا ملع ربق بتاع7هک1 200 ب) 06 وس).گ مظ6006۲50آگیت ۸ م0265 
321-5328)1974(۰ 00۰ و8 ر.۳۱0۳1 ۷۵۵/۰ ۷۱۵۸۱۱۵۲۰ ۰ ۲۱۲ رفعتنتن‌ناطه 

۶ زو هه ونازدانع۱] عطا 10۲ ععوو عل تفصتمتطنه فصتاحمصهع 1۵۲ هعمج اصعت تاه هر ربق رتاع 1۲2۷ 
.(339-344)19860 00۰ ر2 ۷۲6/۲۰ ,تعصصتا ,اممه ,امن .۰ فعتنت ناه 

مامت .فتووآممه امستهناتی ۵۶ 0مطاعصه م۲۵۲6 عمط تم قعوهها عم افصصتصتططانای ریق رتم1۲2۷ 
277-280)1987(۰ 0۳0۰ و3 1۷۲6۰ .م۲ 

413-420)1986(۰ .00 ,24 تاک تام .۱60600 )لهج کم ممتتمع ت0۵ ربق رطع ۲2۲۷ 

رتاک تاوصم .من تصعصعام له اههمه تم ماویری عتقصتتمم ازتاعصم ۸ ربق رتم1۲2۷ 
۳65)1991(۰ 

تناو اقاعاعای ۵۶ فتورلمصه بوتانطانع عمط ۶۵۲ قطاموفتم ۵۶ وعوعده عل7 لقصصتامهطناو "یم رتع1۲2۲۷ 
239-271)1988(۰ 00۰ رك7 .۲۳۱808۵ ۱۷۵0۲۰ ,موه م۱۷ ناهن 

۰ 19 2۷۵۵۸۷۵0۲۱6 .ل ۲ ,و6۲0۲ ۵۴ وتولامصه طقمص ب۶۵ فطاوهم ۵۶ وععده علع7 افصتاوه ۰۵2 ربق رتاع 1۲2۷ 
.(273-279)1989 

,69 7۸۷۲۷۲ ."5ع16اف نازانه۱؟ ععتهوی تم عطاوهتم ۵۶ وععده عل افصتامه 0ص امصتصته ۵ ربق رتم1۲2۷ 
.(214)1989 -212 0۳۰ 

:)۵۱و۲7 را قطموع ما مومت تقصمعی طان من ۷6۵۵۲ 0۲ فصعاوو عتقوظ ر. .6 5602۳66 
164(۰) (۳۱۹51۵0 171-17510 0۴۰ ,19 رانا 

٩۱55120( )1969(‏ جر ) رفن۱۵۲0 قافن ر "قطوجه منتصز۴ ۵۶ ۰16۵۲ رب بظ ,7710۲ 
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سال بیست و دوم» شماره یک ۱۳۸۹ نشریه مهندسی عمران دانشگاه فردوسی مشهد 


1 استفاده از روش‌های نظریه گراف‌ها برای بهبود ... 


عع۷ ناهن رد ععووط ع0۲0 260تامتعصعم امصامن امنهر ومتاهمتصتاوه اقتتمافصتطاصوی ۸ ریق رطع1۵۷ 
.(39-32)1979 0۳0۰ ,20 ,.۹۳۵0۲60 ۷۲۵۵۲۰ ,۸۵۵1 

1980(631-۰) 1۷]2).,16 .22605 ,طاححتع ه ۵۲ وتعوط 6708 مات و ن یط رقافطتعها۲۵ 

بل 51۸۲ ,طححتع و ۵۶ ونعهها ع70 )عمط عطا 4ص ما طاتمعفاه عصنا. اقتجمصرامج ۸ ردنا یگ رط۲0۲0 
358-366)1987(۰ 00۰ و16 رءلا00) 

(2006) ,کل راعکت۹03 :(6) رملز۷۷ ۲۰(ص0ع 2۳۶) عتی‌امصه آهس‌نی آممهتام0" ربق رطع162۷ 

.(1992) , 1,0800۳ ,16۳ ,"کلمطامحظ یتتاه۱۷ مصه امرهتی ۰ :ممتصمطع] ۱۷ امتستتام‌تتتاگ ربق رتم۲۵۷ 
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